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psj ■ tational wave detection. An issue in wavelet analysis is what is the time-frequency 



We present a new method of wavelet packet decomposition to be used in gravi- 



resolution which is best suited to analyze data when in quest of a signal of unknown 
shape, like a burst. In the other wavelet methods currently employed, like LIGO 



WaveBurst, the analysis is performed at some trial resolutions. We propose a de- 
composition which automatically selects at any frequency the best resolution. The 

> , 

0^ , criterion for resolution selection is based on minimization of a function of the data, 

I named entropy in analogy with the information theory. As a qualitative applica- 

tion we show how a multiresolution time-frequency scalogram looks in the case of 
. a sample signal injected over Gaussian noise. For a more quantitative application 

of the method we tested its efficiency as a non-linear filter of simulated data for 

•rH , 

^ ' burst searches, finding that it is able to lower the false alarm rate of the WaveBurst 

algorithm with negligible effects on the efficiency. 

PACS numbers: 04.80.Nn,04.80.-y,07.05.Kf 



I. INTRODUCTION 



It is predicted that gravitational waves are emitted by a variety of cosmological and as- 
trophysical sources and several detectors are now operating, or soon will, to directly observe 
a gravitational signal. Among all possible kind of signals, bursts are generically expected to 
be produced by compact astrophysical sources, like supernova explosions or compact object 
merging, as well as by high energy phenomena, like gamma rays bursts or more exotically 
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by cosmic string cusps and kinks. Not all of these processes are sufficiently well understood 
to predict a specific gravitational waveform. We use here "burst" generically, to designate 
a not better specified transient signals, with short duration, typically less than a second, 
and unknown shape. The identification of bursts, possibly drowned in the detector noise, is 
then an important goal in gravitational wave data analysis. Solid detection algorithms are 
required to generate trigger of events out of the time data stretch coming out of a detector. 

Time- frequency decompositions for analysis of time series have been used since last decade 
in other research field than physics, see e.g. [9}] for an application on pattern recognition 
and [3'] for an application to physiological time series. More recently they have been applied 
to the gravitational wave detector data analysis, see e.g. . In particular the WaveBurst 
algorithm 0, 0] has been employed in the the LIGO burst search, see also 6] for a wavelet 
method based on a different statistics than WaveBurst. 

Here we focus our attention on the issue of selection of the time-frequency resolution. 
Since the signal one is looking for is unknown, it is a priori not clear what will be the reso- 
lution enabling to gather the power of the signal in the smallest possible number of wavelet 
coefficients, and consequently to better single it out from the noise. We then implemented a 
method to automatically recognize from the data themselves what are the best resolutions 
allowing to concentrate the power in the data in the least possible number of wavelet coeffi- 
cients. That is obtained by considering an entropy function, defined so that it is maximum 
when the power is democratically distributed among all the coefficients and minimum when 
all the power is concentrated in only one coefficient. 

The plan of the paper is as follows. In Sec. [Ill we explain how we applied the adaptive 
multiresolution method, showing a case study of a signal injected over Gaussian noise and 
how it is represented in the wavelet decomposition with varying time-frequency resolution. 
In Sec. mil we use this method as a ffiter, showing that in combination with the WaveBurst 
algorithm it can be used to lower the false alarm rate, leaving almost unaffected the detection 
efficiency of the WaveBurst trigger generator. We summarize our conclusion in Sec. llVi 



II. THE WAVELET DECOMPOSITION AND THE ENTROPY CRITERION 

We use a wavelet packet transform to trade a set of data represented by a time series 
with another set defined in the time-frequency domain. 
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The output of a wavelet packet transform can be represented by a binary tree where each 
node, labeled W^, is an orthogonal vector subspace at level j and at layer i. Each of these 
levels will contain exactly the same information of the original time series, so that complete 
reconstruction of the original signal is possible from any level. 

For instance, if the time resolution of the original series, made of 2" samples, is Atg, 
information can be obtained from zero frequency up to a maximal frequency fmax = (2Ato)~^ 
(we restrict here for simplicity to the case in which the number of sample is a power of two, 
but the method is completely general). Each level contains 2"^ coefficients, arranged in 2^ 
layers, each of which has 2"'~-' coefficients. At level j the time resolution is Atj = 2-^Ato 
and the frequency resolution is Afj = fmax- The coefficients in the i^^ layer refer to the 
frequency bin characterized by i2^^ < f / fmax < {i + 1)2"-'. 

Only a subset of these subspaces are needed to completely reconstruct the original time- 
signal: it can be shown (see e.g. 3]) that the set of leaves of every admissible tree completely 
represents the original signal in the wavelets domain, where an admissible tree is a sub-tree 
of the original binary decomposition tree where every node has exactly or 2 children nodes 
(see e.g. fig. [H where the leaves of an admissible tree are marked by red circles). 
Having such a redundancy, how to choose the set of to represent the signal? In ^ 
the are chosen setting j = L, where L is the decomposition level and consequently 
i = 0, ...,2^-^ Here we present a way to choose the based on an entropy function 
E'^{W-), where W- is the set of the coefficients of the base of a node in the decomposition 
binary tree T, and defined as follows 

^W) = -E^iog,^, (1) 

where the Xk G TV/ are the wavelet packet transform coefficients of a Wf node and 

\\x\\'=j:^i (2) 

i 

From the complete decomposition binary tree T of a wavelet packet transform, we select 
the admissible tree Ae C T whose leaves G A have the minimum cost, i.e. where 
E'^i^i ) is minimum for all the children nodes of in the original wavelet packet transform 
binary tree T [entropy criterion Sj, where a different notation is used). 

From among all the admissible binary trees, Ae is the one that represents the signal most 
efficiently. By "efficient" we mean that a signal can be represented by a small number of 
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wavelet packets, that is, the basis for the decomposition is chosen such that the weight of 
the coefficients is concentrated on a small number of wavelet packets and a large number 
of coefficients are close to zero. In fact from eq.([l]) we can see that if in a particular basis 
the decomposition produces all zero coefficients except one (i.e., the signal coincides with a 
wavelet packet transform wave form), then the entropy reaches its minimum value of zero. 
On the other hand, if in some basis the decomposition coefficients are all equally important, 
say Xk = 1/N where is the length of the data, the entropy in this case is maximum, 
log2 A^. Any other decomposition will fall in between these two extreme cases. In general, 
the smaller the entropy the fewer significant coefficients are needed to represent the signal. 

As an illustration we show in figs. [2] and [3] two time-frequency decompositions, or scalo- 
grams. In fig.[3]we report an example of multiresolution scalogram, where we consider a time 
series with a sampling time of 0.2msec, and then performed a wavelet transform up to level 
9 using our adaptive multiresolution algorithm which automatically detects the resolution 
which is better capable to "concentrate" the signal. As a comparison in fig. [2] we show the 
same data in a time- frequency analysis done at decomposition level 9. 




FIG. 1: Schematic example of a binary tree for a Z = 3 wavelet decomposition. The layers in each 
level are explicitly labeled for I < 3. The red circles show an example of a possible choice for 
wavelet packet basis. In this examples it is made by Wq, W2-, W^, W/. 
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FIG. 2: Scalogram in the time-frequency plane. The diagram shows white noise data with an 
injection h{t) = Ae''^^ /^'^^^hm{2iTf), with to = 6.4sec, A = 0.06sec, / = 930Hz and A ~ 0.76. 
The wavelet transform is made on a Symlet basis with maximum decomposition level equal to 9. 

III. TEST OF THE DETECTION ALGORITHM 



Within the context of a wavelet packet transform several detection procedures has been 
developed (see e.g. Q,3 and fl). Here we adopted the WaveBurst method as in [4] and we 
show a possible application of the entropy criterion to improve the method. The WaveBurst 
algorithm is a trigger event generator, it identifies candidate events by looking for power 
excess in a wavelet decomposition and collecting clusters of higher coefficients. We suggest, 
as a possible application of our entropy multiresolution method, that a signal can be filtered 
by looking for excess power not in the wavelet decomposition with some given resolution 
(i.e. by choosing a basis belonging to the same decomposition level), but in the wavelet 
decomposition which makes use of the entropy method, see fig. [H After having zeroed 
the coefficients under a given threshold, one can reconstruct from the surviving wavelet 
coefficients a new, filtered, time series which can be fed to the WaveBurst algorithm. 

The original and the filtered data are both time series and the WaveBurst algorithm can 
be run on both of them. We expect that if a signal is present in the original series, it should 
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FIG. 3: Scalogram in the time- frequency plane as in fig. [2l but from entropy multiresolution 
method. 

be stable under this analysis, or better, that the entropy-mult iresolution method can catch 
all of the signal with a different noise distribution over the wavelets coefficients. So we can 
check the coincidences between the event found by WaveBurst on the original time series 
and on the entropy-filtered one. 

We simulated one hour of white noise data, to which we added 71 injections separated 
by 50 seconds. The injections had the temporal profile 

h{t) = Aexp (^-^) sin(27r/t), 

with / = 930Hz, At = 0.06 sec and A ~ 0.76. 



We first considered how many injections were recovered (efficiency) and how many false 
alarms were found by applying the WaveBurst method to the original data to draw a Receiver 
Operating Characteristic (ROC). We then considered the above mentioned coincidences 
between the events found on the original data and the events found on the data which had 
undergone our adaptive entropy multiresolution filter. The two ROC curve are shown in 
fig. m and they show how the coincidence method enables to lower substantially the false 
alarm rate without affecting the efficiency. 




False Alarm Rate [Hz] 

FIG. 4: ROC curves for WaveBurst method applied on a time stretch of white data with sine 
Gaussian injections. The blue squares refer to the straight WaveBurst algorithm, the red circles 
to the coincidences between the events found by WaveBurst on the original data and on the data 
filtered after the adaptive multiresolution decomposition. On the ROC curves the percentage of 
coefficients which are zeroed varies. 



IV. CONCLUSIONS 



We presented a new method of wavelet packet decomposition to be used in gravitational 
wave detection. This entropy multiresolution method is aimed to better represent in the 
wavelets domain a signal embedded in noise in order to built a better filter in this domain. 
We presented also one of the possible applications of this method, showing how to use it in 
order to have an improved ROC curve by well tested methods, like WaveBurst. In following 
work we hope to present other applications of this entropy multiresolution method. 
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